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Recent experiments by Mertz et al (Phys. Rev. B55, 9160, 1997) demonstrated that the loss of 
superconductivity upon Pr doping is associated with the change of the character of the oxygen holes 
from p a to p n . This experiment sheds new light onto the long-standing problem of T c suppression 
by Pr, and helps to rule out a number of theoretical models, leaving only those which predict such 
a transfer of the O holes. To distinguish between the two models that do predict such an effect, one 
has to access the ratio of the planar and axial character of the holes. We do so in this paper. 



One of the most exciting cases of superconduc- 
tivity suppression in High-T c cuprates is that of 
i?-Ei_ K Pr 2: Ba2Cu307, where RE stands for a rare earth. 
A numbet. of explanations, including magnetic pair 
breaking!!™, disorder on Ba siteB, and hole transfer from 
planes to chainso, have been suggested. All these models 
have problems of various degrees explaining the available 
experimental data. Most notably, recent polarization- 
dependent O Is near-edge X-ray absorption measure- 
ments on single-domain Pr-doped single crystalsB suc- 
ceeded in measuring not only the total amount of O holes 
in CuC>2 planes, but also the relative weight of their pla- 
nar (jp x , p y ) and axial (p z ) components. It is hard to 
overestimate the importance of this experiment, which 
essentially rules out the models that assume the total 
hole concentration in planes to be dependent on the Pr 
doping. 

To the contrary, this result strongly favors the models 
that assume that hole depletion of the superconducting 
pda states occurs because of the hole transfer into the pdir 
states. The first model of this sort is that of Fehrenbacher 
and Ricetj. They assumed that p^ orbitals of the planar 
oxygens neighboring a Pr ion, form a hybride state with 
Pr fz(x 2 -y 2 ) orbital. They noticed that oxygens around 
a rare earth site form a nearly perfect (better that 10%) 
cube, so the lobes of the f z ( x 2_ y 2\ at a rare earth site 
point nearly exactly towards the oxygens. They assumed 
further that such a hybride orbital is strongly localized. 
This assumption followed from the implicitly implied in 
Ref. ^ statement that the corresponding combination of 
the oxygen orbitals forms only when the rare earth site is 
occupied by Pr, but not by Y (where the f z {x 2 -y 2 ) state 
is unoccupied and lies much higher the O p states), in 
other words, that direct hopping between involved O or- 
bitals, without assistance of the Pr f z ( x 2 -y 2 ) orbital, is 
negligible. This seemed to be a plausible idea, because 
in the considered Os cube the 0-0 distance is relatively 
large, do-o ~ a/V% ~ 3 A. Besides, the hopping matrix 
elements between two orbitals pointing to the center of 
the cube is to-o = + ^t pp7T , more ppn than ppa. 

On the other hand, closer look at the band structure 
shows that the 0-0 hopping may be not small even in the 
absence of Pr. For instance, Andersen et alu found sub- 



stantial ppa hopping perpendicular to the planes (that 
is, between p z orbitals), of the order of 0.35 eV. 0-0 
hopping inside the planes is even larger, because it is 
assisted by the diffuse Cu s orbital (see Ref. [?] for the de- 
tails). Thus, it is a priori unclear whether the p^ orbitals 
should form a localized state or rather a band compa- 
rable in width with the superconducting ppoj-.band. To 
answer this question, Liechtenstein and MazinB identified 
among the manifold of occupied O-derived bands the oxy- 
gen band that has the z(x 2 — y 2 ) symmetry, and found it 
to be rather dispersive (band width sa 1.5 eV). Based on 
this finding, they suggested a modification of the origi- 
nal Fehrcnbachcr-Ricc model, in which holes are trans- 
ferred from the pda band not to a localized state, but 
to the band derived from the Fehrenbacher-Rice states 
with substantial overlap between them. Parameters of 
the model were found from LDA+Up r calculations, which 
should be valid at small and moderate Pr doping (at 
large doping, that is, close to the metal- insulator tran- 
sition, one has to take into account Hubbard U at the 
Cu site, which was not done in Ref. |J). This modified 
Fehrenbacher-Rice model was very successful in describ- 
ing T c suppression rates in i?i?i_ x Pr 2: Ba2Cu307, includ- 
ing the dependence of the suppression rate on the rare 
earth. The latter constitutes a serious problem in the 
original modelD , because the only effect of substituting 
on RE by another in that model is "chemical pressure" , 
that is, contraction of the Cu-0 and RE-0 bonds. The 
two contractions shift the energy balance in favor of the 
pda band and the ppn — f z ( x 2 -y 2 ) state, respectivelya. 
However, the latter effect is stronger (cf. canonical scal- 
ing of Harrison, which implies the d~ 4 and d~ 5 scaling 
for pd and pf hopping amplitudes, respectively). Corre- 
spondingly, external pressure enhance T c suppression, the 
most striking manifestation being the recently observed 
suppression of superconductivity by external pressure in 
pure NdBa2Cu306.7t3. To the contrary, bond contrac- 
tions due to host RE substitutions decreases the T c sup- 
pression rate, opposite to what one would expect in the 
localized model of Rcf.0. On the other hand, in the 
modified model of Ref. gthis effect appears naturally as 
a result of hopping between the host RE and the Pr / 
orbitals. 
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However, closer to the metal-insulator transition the 
model of Ref. 0, becomes less and less realistic and it 
was pointed outB that uncritical extention of this model 
to pure PrBa2Cu307 leads to serious problems with ex- 
plaining away possible metallic conductivity in the ppir 
band. Thus, it is highly desirable to address the differ- 
ences between models of Refs. ^and [| by an independent 
experiment. This was one of the goals of Mertz et aO. 
Essentially, they noticed that the p^ holes have different 
character in the two models: According to Ref. ||, they 
are predominantly planar (p x ,y), while according to Ref. 
U the pn orbitals in question point towards neighboring 
Pr, and correspondingly have comparable amount of the 
p x , y and of p z character. 

As mentioned above, around each Pr eight nearest 
neighbor O ions form nearly perfect cube; correspond- 
ingly, p orbitals pointing to the Pr according to the model 
of Ref. ^ form an angle a — arctan(l/\/2) » 35° with 
the xy plane (in reality, the interplanar 0-0 distance is 
slightly larger than intraplanar one, so the angle is ss 36° 
for 02 and f» 37° for 03). Incidentally, in the origi- 
nal papem this angle was mistakenly identified as 45°, 
which dicLnot influence the calculations of Fehrenbacher 
and Ricc,El, but misled the authors of Ref. |^ and affected 
their comparison of the experimental results with the 
Fehrenbacher-Rice model. Furthermore, as it was men- 



tioned in Ref. |8j, in the band model the ppir hole states 
are strictly planar only in the limit when the number of 
ppir holes tends to zero, i.e., when the ppir band touches, 
the Fermi level. The experimental results of Mertz et a£l 
correspond to substantial Pr doping (x = 0.2) and con- 
siderable number of the ppir holes (0.25 holes/cell). In an 
elegant analysis of their experimental data, Mertz et al 
deduce that the angle that unoccupied p^ orbitals form 
with the xy plane is 20°-25°, and conclude that this is 
equally far from the Fehrenbacher-Rice prediction of 45° 
and of the Licchtenstein-Mazin prediction of 0°. However, 
as we mentioned above, the former model predicts 35-37° 
instead, and the latter predicts an angle which has yet to 
be calculated. Below .ive present such a calculation. 

The original paperu was aimed at calculating T c sup- 
pression rate at small Pr concentrations, and utilized a 
tight-binding model where p z orbitals were not included. 
In order to analyze the hole character at finite dopings, 
we have to include both p x ^ y and p z states. We will, how- 
ever, neglect the pfir hopping and leave only the pfa 
one. The hopping amplitude between the f z ( x 2 -y 2 ) an d 

the p XiV orbitals is ^ \f\tpf<? = y§ ^ an d between the 

fz(x 2 - y 2 ) and the p z it is = J\t (Ref. HJ). 

Correspondingly, instead of the tight-binding Hamilto- 
nian given in Ref.H, Eq. (1), we have 
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where the first two orbitals are p z and the remaining part is the p x — p y — f Hamiltonian used in Ref. The same 
notations as there are used, namely S x , y = 2 sin(ak x , y /2), C x . y — 2cos(afc XlJ/ /2). 

In Ref. |^ we were interested in calculating the doping dependence of the number of unoccupied states, that is why 
we had to keep the / orbitals explicitly in the tight-binding Hamiltonian. Now we want to estimate the average weight 
of the p z orbital compared to that of the p XtV orbitals for a given number of holes. Thus, we fold down the / state in 
Eq. ([[]), using Lowdin perturbation theory. This gives 
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where we put e p to zero and denote t — t 2 /3(Ep — ef). 
We diagonalized this Hamiltonian numerically and cal- 
culated the average angle for unoccupied (hole) states 
in the upper p n band as a function of filling (Fig.l) 

(a = arctan ; where n x + n y and n z are cal- 

culated by adding squares of the first two and the last 
two components of the eigenvectors). For the hole count 



0.25, corresponding to the experiment of Ref. |5J, we find 
a = 18°. Fig.l was calculated with r = 0; for fixed hole 
count, the angle a depends on r very little. For instance, 
for r = t we find a to be reduced by about 20%, that 
is, to 15°. The latter number can be taken as the lower 
bound for the band model prediction. In fact, the reason 
that the band model yields an unexpectingly large rota- 
tion angle is that although this angle tends to zero for 
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infinitesimally small doping, it increases very rapidly, as 
the square root of hole number and, consequently, as the 
square root of the Pr concentration, the fact that was not 
appreciated in Ref. |^. 

To conclude, we have calculated the relative charac- 
ter of the planar, p x ,y, and of the axial, p z , hole states 
in the ppn band in Yi_ :E Pr :E Ba2Cu307 as a function of 
hole count in this-, band, using the localized model of 
Fehreabacher-RiceO and the band model of Liechtenstein- 
MazinEi. Comparing our results with the X-ray absorp- 
tion experiment of Mertz et aM, we observe that, in agree- 
ment with the conclusion of the authors of Ref. ||, their 
result falls between the two theoretical predictions. How- 
ever, the two predictions are closer to each other than 
what was supposed in Ref. ^, and correspondingly closer 
to the experiment. Especially close to the experiment 
(10-30%) are predictions of the band model of Ref. ||. 

The last remark we would like to make is that, al- 
though the quantitative difference between the two model 
appears to be not that large in this case, there is a quali- 
tative difference, which can be used in the future experi- 
mental studies to distinguish between them. Namely, the 
band model predicts that rotation angle increases with Pr 
doping, especially at low doping levels, while according 
to the localized model it should stay essentially constant. 
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FIG. 1. Average rotation angle of the depleted 0(2,3) 
Ptt orbitals as a function of the number of holes, n, in the 
ppir ( "Fehrenbacher- Rice" ) band, from the band model of 
Liechtenstein and Mazin (LM). Experimental estimate at 
n — 0.25 is shown by the double-headed arrow, as well 
as (independent of the hole count) rotation angle in the 
localized Fehrenbacher- Rice model (FR). The LM curve 
corresponds to r = in Eq.(^). 
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